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"£ji This paper deals with the reconstruction of both the internal geometry and the 

S values of an inhomogenous acoustic refraction index through a piece-wise constant 

i iterative approximation. We propose to enhance a classical Gauss-Newton scheme 

through the use of geometrical informations gained from a previously developed de- 
fects localization method. Both enhancements are aimed to reducing the number of 
computed parameters. The first one consists in a preliminary selection of relevant 
parameters, while the second one is an adaptive refinement to enhance precision with 
a low number of parameters. Each of them is numerically illustrated. 

1 Introduction 

We consider the inverse medium problem consisting in the reconstruction of an inhomogeneous 
acoustic refraction index from far-field measurements generated through plane waves. The inverse 
medium problem is non-linear and ill-posed, which makes it difficult to solve [11]. Applications 
are, for example and among others, non-destructive structures testing or biomedical imaging 
b [15, 17, 12]. 

Classically, the inverse medium problem is written as a least-squares problem and we will 
look for the index of refraction as a piece-wise constant function. We are thus in the context 
of non-linear and ill-posed parameter identification, which calls for specific methods. Among 
them, for its ease of implementation and efficiency for medium-sized problems, we consider the 
popular Gauss-Newton method applied to a Tykhonov regularized data misfit cost functional [8]. 
However, in its usual form, the Gauss-Newton method treats all parameters alike. Yet, during 
the reconstruction, or even right from the beginning, some parameters can reach an acceptable 



precision while the others still need to be updated to reach a global precision. In the absence of 
some local information, these acceptable parameters are then uselessly updated at each iteration. 
In this setting, we thus explore two uses of geometrical information, gained through defects 
localization, to focus the reconstruction and consequently lighten this inverse medium problem's 
numerical resolution. 

Defects are generally understood as perturbations with respect to some known reference. Con- 
versely, this (known) reference can also be seen as an inexact estimation of the (unknown) per- 
turbed state. It has recently been proved in [10] that the location of the defects in a given 
refraction index could be established from far-field measurements of the unknown state and 
computed through a fast numerical method. We propose here to use this fast local informa- 
tion to reduce the computation effort in the complete refraction index reconstruction process. 
Coupling between shape reconstruction and parameters identification has also been proposed 
by using the Linear Sampling method [5] in [4]. However, the Factorization method approach 
used [10] provides a more straightforward formulation as an equivalency that is defined at each 
point. 

In a first time, we consider the case where a known index has been locally modified. This 
could happen, for instance, from a local deterioration or a partially incorrect estimation of the 
actual index. In this case, a preliminary defects localization provides a geometrical information 
that is which parameters have to be reconstructed. The reconstruction can then be performed 
straightforwardly on a reduced computational domain. In a second time, we investigate adap- 
tive refinement. Here, defects localization is used to exhibit inaccurate regions in the current 
reconstruction. This local information allows us to add base functions in selected locations and 
resume the reconstruction to get a better precision while restraining the number of unknown 
parameters. 

This paper is structured as follows: At first, section 2 specifies the mathematical setting. The 
direct problem is introduced in section 2.1 and followed in section 2.2 by the description of the 
inverse medium problem we are interested in. The numerical method on which we will build 
our enhancements is then described in section 2.3. Afterwards, the defects localization and its 
applications are presented in section 3. We show how to reduce the reconstruction domain in 
section 3.1 and the adaptive refinement process is detailed in section 3.2. Finally, we numerically 
illustrate the coupling of both applications, and also on a non-trivial example, in section 4 and 
conclude. 

2 Presentation of the problem 

We here introduce the inverse medium problem, followed by its numerical treatment and preceded 
by the direct problem. 

2.1 The direct problem 

If we consider time-harmonic acoustic waves with a fixed wave number k, the spatial part of 
the wave equation is modeled by the Helmholtz equation [6]. Inhomogcneous media are then 
represented by an acoustic refraction index, denoted by n <E L°°(M. d ), and so, the total field, 
denoted by u n € Lf oc (R d ), is assumed to satisfy 

Au n + k 2 n(x)u n = 0, x e R d , (1) 

where d is the problem's dimension (d = 2 or 3). We consider compactly supported inhomo- 
geneities and denote by D the support of n(x) — 1. Also, we denote an incoming wave satisfying (1) 
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Figure 1: General setting and notations. 

with n = 1 by u l £ Lf oc (R d ). The total field is then the sum of this incoming wave and the wave 
scattered by the inhomogeneous medium, denoted by u s £ Lf oc (R d ): 

u n :=u s +u\ (2) 

where the scattered wave is assumed to satisfy the Sommerfeld radiation condition 

d r u s =iku s + o(\x\~^y (3) 

Then, the linear system (l)-(3) defines u n uniquely from it 4 , and it is known to be invertible in 
L 2 (D) [G, Chap. 8]. 

Besides, the outgoing part of a wave has an asymptotic behavior called the far field pattern, 
denoted by u^° £ C°°(T m ) and given by the Atkinson expansion [19] 

u n {x):=u\x)+ 1 ^-^ T u™(x)+o(\x\- ± ^\, x:=^-£T m , (4) 

ja-jT- V J \x\ 

where T m denotes the set of measurement directions as a subset of the unit sphere S d ^ 1 (see 
figure 1) and where 7 depends only on the dimension and is defined by 

f <^JL if d = 2, 

Furthermore, for practical reasons, we will mainly consider scattered waves having a plane- 
wave source. These plane- waves are defined by 

?/(#, x) := exp(ik9 ■ x), 

where 9 is a unitary vector in the set of incidence directions, denoted by r e as depicted on 
figure 1. We then denote the total field with a plane- wave source of incoming direction 9 by 

u n (9,x), §£T e ,x£R d . 

Lastly, the corresponding far-field pattern in the measurement direction x £ T m will be de- 
noted by 

u™0,x), 9£T e ,x£T m . 

2.2 The inverse medium problem 

We are interested in the reconstruction, from far-field data, of an (unknown) inhomogeneous 
refraction index that will be denoted throughout this paper by n* £ L°°(D). All considered 



indices are implicitly extended by 1 outside D. We thus define the index-to- far- field mapping 
T:L 00 (D)-^C 00 {T e xT m )by 

.F(n):=<°- 

Besides, data are generally perturbed by noise or measurement errors and we only have access 
to u e satisfying 

||u £ -<i|le«>(r e xr m ) see. (5) 

The usual approach to this problem is then to minimize the non-linear data misfit. However, 
most methods to solve this problem are set in Hilbert spaces, so we will have to consider T as 
an application from L 2 (D) into L 2 (T e x T m ) and thus define the data misfit by 

J(n) := ci HJ'(n) - u e ||^ 2(rcXrm) , 

where C\ denotes a normalization constant (e.g. c\ := |u £ || L2 fp xr \). 

Even so, this problem is not continuous, as is shown by the following proposition. As such, 
even a small perturbation e can lead to a minimizer very far from n*. 

Proposition 1. The non-linear problem consisting in "finding n e minimizing J" is ill-posed in 
the sense of Hadamard. 

Proof. The application T is compact and thus cannot have a continuous inverse. Indeed, it 
has been shown that the total field u n is bounded with respect to the L°° (D)-norm of n [16, 
Proposition 2.1.14]. As a consequence, the same holds for the mapping n H> (n — l)u n . Moreover, 
the asymptotic behavior of the Lippmann-Schwinger equation yields the following relation [6, 
Chap. 8.4]: 

J r (n)(0,x)= [ e- ik *- z k 2 (n{z)-l)u n 0,z), 9eT e ,xeT m . (6) 

JzeD 

Hence, the non-linear operator T is the combination of a linear compact operator with a contin- 
uous mapping and is thus compact itself. □ 

2.3 Iterative approximation by a piece-wise constant index 

As stated in the introduction, we try to recover the unknown index n* with help of piece-wise 
constant functions. Hence, the indices will numerically be represented by N complex parameters 
(j/»)»=i...jv associated to the same number of zones (2i)i=i...jv, so n(x) — J2i=i n Vi^-Zi(x). Each 
zone is thus a set of connected elements in the underlying mesh used to solve the Hclmholtz 
equation. Moreover, to avoid any possibility of inverse crime, the reconstruction will be led on a 
different mesh than the one used to generate the data u e . On this behalf, the respective supports 
of the basis functions associated to the unknown parameters will intersect the discontinuities 
of n*. It is thus strictly impossible to reconstruct n* exactly. This is illustrated on Figure 2, 
showing two 2D meshes that will be used our numerical examples. 

In this setting, we use the well-known Gauss-Newton method applied to the cost function J 
with a standard Tikhonov regularization term [13]: 

J(n) := ci ||J"H - u e ||^ 2(rcXrm) + c 2 ||n - n \\ 2 L2{D) , 

where c 2 > is a regularization parameter and no € L 2 (D) is an initial guess. Assumptions on 
no and c 2 for this method to converge are discussed in [1, 3]. The sought-for index n* is hence 
approximated by a sequence (n p ) pe N of indices, defined iteratively through Algorithm 1. 
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Figure 2: Test case geometry 



Input: n £ L 2 (D) 

1 p<-0; 

2 repeat 

Compute rip+i by solving the linear system 

[DFfapYDTinp) + —id)(n p+1 - n ) = 



-DT(n p )*(T(n p ) - u 6 - DT(n p )(n p - n )), 



where id stands for the identity matrix; 
p<-p+l; 

5 until ||TOp-np_i|| 2 /(l + ||n p _i|| 2 ) < e; 
Output: n PEnd 



Algorithm 1: The Gauss-Newton method for J 



Moreover, we consider small enough amounts of data and unknown parameters to be able to 
compute the cost function's gradient in a reasonably fast time on a standard desktop computer. 
This gradient has the following simple formulation. 

Lemma 2.1. The application J- is twice differentiable. Moreover, the differential DJ- evaluated 
at n £ L°°(D) and applied to the direction dn £ L°°(D) has the following integral representation 

DJ-(n) dn : (9,x) H> / k 2 u n (—x,z)u n {9,z)dn(z)dz, 9£T e ,x£T m . (7) 

JzeD 

Proof. Expansion (4) shows that u™(6, •) is linear with respect to the scattered field (u n — 
v?){6, •). Furthermore, It has been shown in [16, Proposition 4.3.1] that the scattered field is 
twice differentiable with respect to n and that the differential of the index-to-scattered-field 
mapping evaluated at n £ L°°(D), applied to dn £ L°°(D), is the function v s (9,-) £ Lf oc (M. d ) 
satisfying 

(A x + k 2 n(x))v s (9,x) = -k 2 u n {9,x)dn{x) 1 x£R d , (8) 

and the Sommerfeld radiation condition. Note that, contrary to n, the direction dn is extended 
by outside D. Thus, T is twice differentiable and its differential is defined on C°°(T e x T m ) by 
DT(n)dn = v°°. 

Now, note <f> n (z,x) the Green function to the Helmholtz equation (1). Multiplying (8) by 
< & n (z, x), integrating over D and using Green's formula then yields 

v s (9,x) = / k 2 $ n (z,x)u n (9,z)dn(z)dz 1 x £ R d . 
JzeD 

The asymptotic behavior is then given by 

v°°0,x)= f k 2 ^(z,x)u n (9,z)dn(z)dz, xGS*- 1 . 
JzeD 

Finally, representation (7) is obtained by applying the mixed reciprocity principle: $^°(z,a:) = 
u n {— x, z) (see [14, equation (3.66)]). □ 

Numerical example 

To illustrate our reconstruction schemes in M 2 , we consider a disc D of radius 1 centered on the 
origin. The embedded perturbation Q is then another disc, this time centered on (0.3, 0.3) and of 
radius 0.3, as depicted on Figure 2a. The (perturbed) index we are looking for is set to n* := 1.3 
in D \ CI and n* := 1.6 in CI whereas the initial guess, corresponding to the last known state, is 
uq := 1.3 in D. 

The Gaus-Newton method is performed with the regularization parameter set to C2 := 10 -2 . 
Also, required Helmholtz equation's solutions are computed via P l finite elements and carte- 
sian Perfectly Matched Layers (PML) [7]. Lastly, the corresponding far-fields are evaluated 
through representation 6. For all these examples, the wave number is set to k — 5 and the 
incoming/measurement direction's angles are equally spaced in [0,2-7r]. 

An example can be seen on Figure 3 with the reconstruction mesh's 2672 triangles (see Fig- 
ure 2b) divided into N = 10, N = 27, N = 75 and N — 2672 zones. For comparison purpose, 
we also list in Table 1 the exact final relative error 

e PE„d := ll n PE„d — n \\l 2 (D) I H n Hl 2 (,D) 

obtained in several configurations. Besides, we see in this table that the relative error is in the 
order of 10~ 3 , so we set the stopping criterion to e = 10 -4 for all our reconstructions. 



0.5 



-0.5 




10" 




-0.5 0.5 

(a) TV = 10 parameters 




-0.5 0.5 

(b) N = 27 parameters 




-1 -0.5 0.5 

(c) N = 75 parameters 



0.5 0.5 

(d) N = 2672 parameters 




T 



^4 



6.5 
6 

5.5 
5 

4.5 
4 

3.5 
4 



-i-N=10 
-x-N=27 

A N=75 

V N=2672h- 



V 



'V 

3 

(f) prt 100e p := 100j|n p -n*||/||n*| 



(c) p h^ J(n v ) := ci \\T(n p ) - u £ | 

Figure 3: Gauss-Newton reconstruction with 30 x 30 data and 2% noise 
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Table 1: Gauss Newton reconstruction 

3 Enhancements of the Gauss-Newton method via defects 
localization 

In the presented piece-wise constant iterative approximation, the possible precision is directly 
linked to the number of parameters N which, in turn, is linked to the computational effort. In 
the lack of some geometrical information, all constants are equally treated and updated at each 
iteration. In some cases however, this can generate more effort than is really needed and we 
address two of them. 

1. For the first case, we place ourselves in the context of a localized perturbation of an 
otherwise known initial state Uq. So, one can state that most of the index has not changed 
and should not be reconstructed. 

2. For the second case, we are concerned in focusing on the most inexact constants during the 
reconstruction. Indeed, to obtain a precise identification, the reconstruction mesh has to 
be refined in the zones intersected by the discontinuities of n*. However, if n* is constant 
in large areas, refining the reconstruction mesh everywhere only raises the computational 
effort for a relatively small precision increment. 

To address these aspects of the reconstruction, the useful information in both cases would 
thus be the localization of the exact (enough) constants. Of course, to enhance the complete 
reconstruction, access to this specific information should be fast. To this end, it has been shown 
that there exist a defects localization function recalled in the following theorem. 



Theorem 3.1. 

operator W := 



[10, Theorem 6.1] Assume T r< 



r e = s 



d-i 



Then, define a measurement 



Id + 2ik\j\ F n ) (F n ± -F n ), where F n : ^{S^ 1 ) -> L 2 {S d - 1 ) denotes the 



classical far-field operator, defined by F n g(x) := (<?, u^{-,x)) L 2t S d-i\. Next, we define the posi- 
tive self-adjoint operator W# by W# := \W + W*\ + \W — W*\, where the notation |-| applied to 
an operator L stands for \L\ :— (L*L)z . Lastly, assume that n and n* are real valued and that 
either (n — n*) or (n* — n) is locally bounded from below on CI :— support (n — n*). 
Then, for each z € M. d , we have 



n(z) 7^ n*(z) 



S{n,n*}(z) ■= Yl 



\(u n (-,z) 1 1p j ) L 2 iS d-l ) 



>o, 



where (<7j, ipj) is an eigensystem ofW#. 

Remark 1. As it is, this theorem requires full bi-static data and real-valued indices. However, 
we also recall the conjecture, stated in [10, Remark 6.2], that the same result could be obtained 
by defining the localization function S with (<Tj, ipj) a right-singular system of (F n * — F n ). The 
main benefit is the possibility to consider T m ^ T e ^ S d ~ l and complex valued indices. 

Furthermore, numerical examples in [10] seem to indicate that this localization is effective 
for defects bigger than (approximately) one sixth of the wavelength. Besides, in order to get 
satisfactory results in the successive resolutions of the Helmholtz equation, we have set the 
average reconstruction mesh element's size to be about one twentieth of the wavelength. Thus, 
we will only consider defects that cover at least four connected mesh elements. 

Finally, the examples shown in [10] exhibit that defects can be located even when the surround- 
ing background is not precisely known. Practically, low amplitude inaccuracies with respect to 
the exact index do not seem to interfere with the localization of the contrasting defects. Thus, 
geometrical information gained through the defects localization presented here is expected to 
focus on the most "defective" zones. 

3.1 Selective reconstruction 

We here consider the case where the initial guess n is exact, except for some perturbation 
whose support will be denoted by fl. Thus, we propose to perform a preliminary selection of 
the parameters that do need to be reconstructed. The whole index n* is then reconstructed by 
updating those parameters only. This should generally lead to a reconstruction using significantly 
less parameters, summarized in Algorithm 2. 

Input: n £ L 2 (D) 
i Si «- max Zi 5{„ 0in *}(x); 

2 Q7- ■£- the set of zones for which Si > Tmax(^); 

3 n PEnd 4— Algorithm l(no|n T ) (&U indices are extended by no outside fir - ); 
Output: n PEnd 

Algorithm 2: Selective reconstruction 



Numerical example 

In the framework of section 2.3, we here use one parameter for each reconstruction mesh element. 
Figure 4 shows which parameters are selected with three threshold values T = 10%, T = 20% 
and T = 30%. We can see that a threshold of T = 10% selects slightly too much mesh elements 
and thus provides a satisfactory reconstruction with only 316 parameters. This is significantly 
less than the initial 2672 and thus, reaches the set goal. More precisely, we can see in Table 2 
that the precision is superior to what was obtained through a full Gauss-Newton reconstruction. 
This was to be expected since, in the present case, all parameters outside the perturbation stay 
fixed to the exact value, where they can be misevaluated in the full reconstruction. 

However, a threshold of T = 20% seems too high, as the 178 selected parameters do not com- 
pletely cover the perturbation's support, resulting in a slightly flawed reconstruction. This brings 
up the problem, arising with every adjustable parameter, of how to select a correct threshold, 
depending at least on the measurements noise and the amount of data. Unfortunately, for the 
moment, we do not have a reliable indicator to tell if the selected threshold was acceptable. 
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Figure 4: Selective reconstruction with 30 x 30 data and 2% noise 
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Table 2: Selective reconstruction 
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3.2 Adaptive refinement 

As stated in section 2.3, since we use a reconstruction mesh that is different from the one used to 
generate the data, these zones will most likely intersect discontinuities of n*. Even more so with a 
low number of zones N. Thus, we propose to iteratively refine the reconstruction mesh with help 
of the previously introduced defects localization, in order to provide a satisfying approximation 
of the unknown index with a small number of parameters. The refinement's outline is presented 
in Algorithm 3. 

Input: no G L 2 (D) 

2 repeat 

3 Si <- max Zi S {rlp ^ } (x); 

4 I 4— {i such that Zi contains more than 16 mesh elements}; 

5 ispiit <- i such that S iSplit = ma,x ieI S t ; 

6 Update the set of zones by splitting Zi g ut into four sub-zones; 

7 Update the set of parameters accordingly by duplicating ?7j g Ht three times; 

8 N ^ A + 3; 

9 "p+p E „d <~ Algorithm \{n p ); 

10 P^P + PEnd] 

n until N > N max or each Zi contains less than 16 mesh elements] 
Output: n PEnd 

Algorithm 3: Adaptive refinement 

The number of 16 mesh elements is taken so that, after the splitting, each zone still has more 
than four elements, which is the lower limit for defects to be relevant, as specified in Remark 1. 

Numerical example 

In the same conditions as before, we illustrate our adaptive refinement on Figure 5. Figures 5a- 
5h depict alternately steps 3 and 9 of Algorithm 3 and show how the reconstruction focuses on 
the contrasting perturbation's support. Figure 5i represents the values of rc.42, which is obtained 
with N = 46 parameters placed during 15 successive adaptive refinements. Also, the relative 
error e p , obtained in step 10 of the algorithm, is plotted on Figure 5j as a function of p. 

Comparing with the relative errors listed in Table 1 , or on Figure 3f , we can see that the guided 
adaptive refinement achieves better results than randomly placed basis functions. This example 
shows results even comparable to the full reconstruction with 20 times more parameters. It thus 
pictures a satisfactory reconstruction with a limited number of well-placed parameters. 

These results present sensitivities to noise or data amount that go in the same lines as the 
previous examples, as can be seen in Table 3. 
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Table 3: Adaptive refinement 
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Figure 5: Adaptive refinement with 30 x 30 data and 2% noise 
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4 Combining both strategies 

The selective reconstruction is presented in section 3.1 as a preliminary step to the reconstruction. 
It can thus be followed by the adaptive refinement described in section 3.2. So, we further 
illustrate selective reconstruction and adaptive refinement by chaining them, as described in 
Algorithm 4. 

Input: n G L 2 (D) 

1 Si <- maxx. 5 { „ 0! „* } (a;); 

2 Vtf «— the set of zones for which Si > Tmax(Si); 



3 n n 



<— Algorithm 3(no|o r ) (all indices are extended by no outside fi-r); 



Output: 



Algorithm 4: Selective reconstruction followed by adaptive refinement 



Numerical example 1 

In line with section 3, we illustrate Algorithm 4 with the selection thresholds T = 10%, T = 20% 
and T = 30%. The respective selected mesh elements can be seen on Figures 4a-4c. Figures 6a-6c 
show the reconstructions after 2, 3 and 17 adaptive refinement loops with a threshold T = 10%. 
As expected, since we already saw that the whole perturbation is captured, the result is a good 
reconstruction. In fact, the exact values listed in Table 4 show that this reconstruction achieves a 
precision in par with the results obtained through the initial selective reconstruction with about 
10 times more basis functions. 

Moreover, again similarly to the examples presented in section 3.1, T ^ 20% results in a too 
small selection, leading to a flawed reconstruction. It can be seen on Figures 6d-6i that the 
reconstruction tends to a crown shape. The corresponding relative error values are presented on 
Figure 6j and detailed in Table 4. 
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Table 4: Selective reconstruction chained with iterative refinement 



Numerical example 2 

As a last example, we now consider a less trivial and complex valued unknown index n* , shown 
on Figure 7. Besides, we also make this reconstruction more challenging by reducing the mea- 
surements aperture. The incoming directions are still taken in [0, 2ir], but there will be five 
less, and the measurement directions are now taken in [0, |7r]. In this situation, the localization 
function presented in Theorem 3.1 can not be defined. So, we consider the technical modification 
recalled in Remark 1 that is conjectured to cover this case. 
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Figure 6: Selective reconstruction chained with adaptive refinement, 30 x 30 data and 2% noise 
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Figure 7: Exact index n* 



Furthermore, we assume that n* was known previously to the central perturbation, which 
translates in the initial guess no depicted on Figure 8. The reference reconstructions obtained 
with Algorithm 1 are presented in Table 5. 
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Figure 8: Initial guess no 



To remain in the previously defined context, we present the results of Algorithm 4 applied to 
this new geometry with the same thresholds T = 10%, 20% and 30%. Figure 9 then depicts for 
each threshold the corresponding selected parameters and the resulting reconstruction. In this 
case, T = 20% now seems to be the best threshold value and T — 30% is still too high. This is 
confirmed on Figure 9j, where we can see that, even though T = 10% allows to reach a satisfying 
precision, it requires much more refinements to do so than with T = 20%. Also, more specific 
results for different simulated configurations are listed in Table 6. 

The results obtained in section 3 are thus backed up by this example, exhibiting reconstructions 
comparable in precision to the full Gauss-Newton reconstruction, but with a much lower number 
of parameters. 

5 Conclusion 

We have used a previously developed defects localization method to propose two ways of reducing 
the number of parameters in the reconstruction on an unknown refraction index. The first method 
is set in the context of defects identification and uses their localization to reconstruct only the 
useful parameters of the whole index. The second method is an adaptive refinement, based on 
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Tabic 5: Gauss Newton reconstruction for the more elaborate example 











10 


X 15 data 










25 


x 30 data 










55 


X 60 data 






r 


e 


N 


PEnd 


e PEnd 


Si 


s 2 


J 


VJ 


N 


PEnd 


C PEnd 


Si 


s 2 


J 


VJ 


N 


PEnd 


e PEnd 


Si 


S 2 


J 


VJ 




5% 


76 


95 


12.6% 


0.89 


0.10 


le-3 


3e-4 


76 


91 


14.9% 


0.70 


0.07 


3e-3 


3e-4 


76 


86 


10.2% 


0.70 


0.04 


5e-3 


3e-4 


0% 


2% 


76 


90 


13.4% 


0.87 


0.06 


5e-4 


3c-4 


76 


84 


9.6% 


0.73 


0.02 


4e-4 


3e-4 


76 


>4 


9.1% 


0.72 


0.02 


8e-4 


3e-4 




1% 


76 


93 


11.9% 


0.87 


0.03 


2e-4 


3c-4 


76 


85 


8.1% 


0.76 


0.01 


le-4 


3e-4 


76 


84 


8.1% 


0.70 


0.01 


2e-4 


3e-4 




5% 


73 


82 


10.0% 


0.89 


0.84 


le-3 


4c-4 


75 


76 


6.3% 


0.70 


0.71 


2e-3 


2c-4 


52 


58 


5.7% 


0.70 


0.61 


5e-3 


2o-4 


10% 


2% 


57 


m 


5.9% 


0.87 


1.19 


4e-4 


2c-4 


48 


50 


5.6% 


0.73 


0.70 


6e-4 


2e-4 


48 


50 


5.5% 


0.72 


0.80 


le-3 


2e-4 




1% 


51 


52 


5.6% 


0.87 


1.32 


2e-4 


2c-4 


48 


47 


5.5% 


0.76 


0.79 


3e-4 


2c-4 


45 


43 


5.2% 


0.70 


0.76 


4e-4 


2e-4 




5% 


51 


53 


5.77c 


0.89 


1.36 


2c-3 


2c-4 


42 


42 


5.5% 


0.70 


0.92 


3e-3 


2c-4 


45 


42 


5.7% 


0.70 


0.75 


5c-3 


2c-4 


20% 


2% 


45 


42 


5.4% 


0.87 


1.44 


5e-4 


2c-4 


34 


31 


5.7% 


0.73 


0.97 


7e-4 


5e-4 


27 


27 


6.3% 


0.72 


0.84 


le-3 


6e-4 




1% 


40 


35 


5.3% 


0.87 


1.48 


3e-4 


3c-4 


24 


24 


6.2% 


0.76 


0.90 


5e-4 


6e-4 


13 


17 


7.8% 


0.70 


0.93 


8e-4 


9e-4 




5% 


37 


35 


5.7% 


0.89 


1.61 


2c-3 


4c-4 


28 


26 


6.4% 


0.70 


1.02 


3e-3 


7c-4 


13 


18 


8.5% 


0.70 


0.88 


6e-3 


le-3 


30% 


2% 


25 


25 


6.5% 


0.87 


1.58 


7e-4 


7c-4 


13 


18 


8.6% 


0.73 


1.08 


le-3 


le-3 


13 


18 


9.2% 


0.72 


0.94 


2e-3 


9e-4 




1% 


16 


19 


7.8% 


0.87 


1.61 


7e-4 


8c-4 


13 


18 


8.9% 


0.76 


1.01 


le-3 


9c-4 


13 


19 


10.3% 


0.70 


1.00 


2e-3 


le-3 



Table 6: Selective reconstruction chained with adaptive refinement for the more elaborate ex- 
ample 
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Figure 9: Selective reconstruction chained with adaptive refinement for the more elaborate ex- 
ample with 30 x 30 data and 2% noise 
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defects localization to iteratively reconstruct a better approximation with a contained number 
of parameters. 

Both methods have shown good numerical results with the presented Gauss-Newton method. 
The reconstruction could however be further enhanced by two automations: some automatic 
thresholding for the defects localization function and some automatic selection of the regulariza- 
tion parameter. The second issue has been reviewed for example in [9, 2] and is claimed to be 
less critical when using a so-called Multiplicative Regularization described in [18]. However, to 
this point, we have not been able to further enhance our results with these techniques. 
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